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摘要 

Toader 平均 与 椭圆 周 长 有 着 相当 密切 的 联系 ， 但 是 数学 家 已 经 证 明了 Toader 平均 不 能 
被 初等 函数 所 表示 ， 只 能 用 初等 函数 来 副 求 上 下 界 . 为 研究 Toader 平均 的 近似 值 ， 本 文 从 
椭圆 周 长 公 式 入 手 ， 引 入 a 、5 的 平方 平均 、 算 术 平 均 、 几 何平 均 、 调 和 平均 ， 以 及 椭圆 离 
心率 上 这 五 个 变量 ， 然 后 用 泰勒 公式 展开 求 得 4 个 平均 和 椭圆 周 长 L 中 龙 各 次 方 项 系数 的 
值 ， 并 作 进 一 步 比较 ， 据 此 推导 椭圆 周 长 的 上 、 下 界 ， 确 定 它 的 近似 公式 ， 从 而 推导 Toader 
平均 的 近似 值 ， 进 而 得 到 一 系列 关于 Toader 平均 的 不 等 式 和 椭圆 周 长 不 等 式 . 

创新 之 处 : 利用 震级 数 展开 和 二 分 法 寻找 Toader 平均 的 最 佳 上 、 下 界 . 

闪光 之 处 : 利用 Toader 平均 与 椭圆 周 长 之 间 的 关系 得 到 了 求 椭 圆周 长 的 较 好 的 近似 公 


式 


关键 词 ，Toader 平均 、 椭 圆周 长 近似 公式 、 平 均值 不 等 式 、 泰 勒 公式 、 项 明达 公式 
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The Toader mean and the approximate 
formula of the elliptical circumference 


Abstract 


Toader mean 1s closely related to the circumference of oval, 
but it has been proved by the mathematicians that Toader mean 
can only be approached by but not be expressed by elementary 
functions. In order to study the approximation of Toader mean, 
this article start from the formula of elliptical circumference, 
introducing five variables: square mean, arithmetic mean, 
geometric mean, harmonic mean and the elliptical eccentricity, 
k, then use Taylor formula to unfold them, obtaining their values 
of the k-th power item's coefficient. Depending on the further 
comparison, it can deduce the upper bound and the lower bound 
of ellipse and determine its approximate formula, so that obtain 
a serious of inequalities about the Toader mean and about the 
elliptical circumference. 


Innovations: Use the expansion of the power series and 
dichotomy to discover the best upper bound and lower bound of 
Toader mean. 

Shines: By the relationship of Toader mean and the 
elliptical circumference, obtain an approximate formula of the 
elliptical circumference with great accuracy. 


Key words: Toader mean, approximate formula of ellipse, 
mean inequality, Taylor formula, Xiang Mingda formula 


Toader 平均 与 椭圆 周 长 近 似 公式 


1. 引言 
对 任意 的 re[0,1]， 本 文 总 记 : r = J1- 2 
对 实 参 数 a,b,c(c #0,-1,-2,:- -) , Gauss 超 几何 函数 定义 为 口 : 


Е(а,Б;с;х) =, F (a,b,c; x) = У «М, х”,хе(-1,1), 
fay (c),n! 


th, (а), = аж0), 


л ШІ 
elr) -ZrO Sif (1— 2 sin2 x)? dx,é(r) =e(r) 


(0<r<1) 称 为 第 二 类 完全 椭圆 积分 . 
作者 在 文 [2] 中 介绍 了 Toader 平均 : 


аав. 
т а 


T(a,b)== [° Va" cos 0+b' sin? 0 40 = aT 1 (2,828: 
л а 


а, a=b. 


其 中 ，& =e(r)=[?v1-7 sin? 0 d6,r e [0.1]. 
文 [3] 中 介绍 了 Holder 平均 : 
1 
a’ +b? \p 0 
M,(a,b)=4\ 2 J >?” 
Vab, p=0 
特别 地 ， 当 万 = 3,2,1,0, 一 1 时 分 别 得 到 : 
a+b 
立方 根 平均 : MENG ONES еті 
2 b? 
PHRA: Sla, b) =, 2, 


2 
算术 平均 ， A(a,b) = = : 


几何 平均 : С(а,Ь) = Jab , 
WAVE. H(a,b) =o, 


a+b 
此 外 ， 还 有 


形 心平 均 ，C(a,b) = 


2(a° +ab+b’) 
3(a+b) 
2 2 
反 调和 平均 ，C(a,b) = tP. 
a+b 


2013 年 ， 在 文 [入 中 ， 作 者 找到 了 下 列 双边 不 等 式 
a,C(a,b) +(1—a@,) A(a,b) < T(a,b) < B,C(a,b) +(1— B,) А(а,5). 


Td “4n=1,2,---IN, (а),-а(а-1)а-2)“(а-п-1),аж0). 特别 地 ， 


C1) 


(2) 


(3) 


满足 


成 立 的 充 要 条 件 是 55, ЕС 
T 
a, 5 а). <T(a b) < —— pb .0- (1-8) N 
A(a, 页 = C(a,b) ` Aa, ТЫ C(a,b) 


成 立 的 充 要 条 件 是 co <л-3, B, 21. 


3 12 ЕСЕК 
(3) Ra,=—, М----3, ЖЕТ: 
4 m 


2 
5а? +6ab+5b° <T(a,b)< 2 = aLi by l 
8(a +b) a+b 


(4) Җа, =л-3, РИ 得 最 佳 不 等 ok: 


ла? +2(6—z)ab+ zb? s ЖЕР 13a? +22ab+13b° 
2(a+b Xa °+ab+b2) T(a,b) 8(a+b)(a* +ab+b’) 
8(a+b)(a° +ab+ b°) 
13а? +22ab+13b2 
束 文 辉 在 文 [5]， 得 到 了 3) 式 的 4 个 等 价 结论 
双边 不 等 式 
a M ..,(a,b)+ H(a,b) <Т(а,Б) < P M, (а,Б) + H(a,b) 


<T(a,b) < 
л 


2 
成 立 当 人 p>—. 
8 A 


双边 不 等 式 
GMa; оа 5)<Т(а,5)< В.М... (а,Б)+ A(a,b) 
ач ча < А 

m 2 


双边 不 等 式 
M. (а,б) + C(a,b) <T(a,b)< B,M cu (a,b)+C(a,b) 


4 3 2 

зу НИХ а, < т bB 24-1. 
T 

双边 不 等 式 

М. (а,Ь) + C(a,b) < T (a,b) < BM (a,b) + C(a,b) 

1 

立 当日 704” <---, SS =. 

成 立 当 且 仅 当 @&， F В, 


我 们 的 主要 结果 有 
sa 各 种 平均 定义 如 上 ， 则 


24128834, b) -579H (a,b) -61G(a,b) -195S(a,b)] < T(a,b) 


ДЫ 
<s> [43.A(a,b) 7H (a,b) -3G(a,b) - S(a,b) 
(12) 式 的 下 界 比 (5) 式 的 下 界 强 ， 上 界 则 不 分 强 弱 . 


a b 


2(a+b)(a* +ab+b’) 
a’ +2(6-л)аЬ+лЬ?. 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


2 2 
由 椭圆 图 形 的 对 称 性 与 Toader 平均 的 定义 可 知 , HL + 2 = Ка > b> 0) WAKL 
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ie [ve cos’ 0 +b° sin? 0 40 = 4[2 Va’ cos’ 0 + b° sin’ 0 40 
= 2n-=[? Va’ cos’ 0+? sin? Ө 10 =27T (a,b). 由 此 可 得 
т 


定理 MEAK L AU MEAR 
TF [2883 A(a, b) -579H (a,b) — 61G(a,b) -195S(a,b)] 


(13) 
<(=)L< (“1143 46а,») ~7H (a,b) -3G(a,b) — S(a,b)] 
2, Be 
为 了 证 明定 理 1. 2, е 
引 理 1 车 4a 2b>0,r=2ef0,1]， 则 
a 
T(a,b)== Va соз 0+2 sin” 40 =“ |° Vr sin? 0 40 
NA 
(2п- Dig pe 
=afl— 
i 2 (2n)!! ЕТЕ г 
1, 3, 5 ¢ 15, 441, 4851 14157 2760615 
SU ааа 8 Рет —— таи 2% p16 — zj 
(14) 
证 明 : 利用 级 数 公式 
es 1 
а уар) с 
2 2 2 Де. querer О, 
Ла > a x 1+2 1) “Omi Хе), 
Ae -pD 
及 [2sin” tar = eo 借助 逐 项 积分 法 可 求 出 T(C, 态 的 无 穷 级 数 表 达 式 ; 
п).. 


Т (a,b) = [2.122 сов? ө +6? зіп? Ө 40 
Л 0 
= 24 |; Vi-r’ sin? 0 d0,r e [0,1] 
z 


= pi» С Jy L 1 or (-r? sin? A)" dt 


= 24 2р Deas DE r” sin?” tt 


(2m)!! 
(2п- DEY р? 
-41- 
t > (2n)!! ro 
1 3 5 175, 441 4851 14157 2760615 
=а1-5” - Sa? oat 7 ace a r? 546 z 26.2 у; 


引 理 2. 若 椭圆 的 标准 方程 是 : ts 4 号 =1a>b>0， 工 为 椭圆 的 周 长 ， 则 有 
a 
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L= 4af? 1-е? sin’ tdt =2xa{l— De De 


(2n)! Беті 
1, 3.4 5, 175, 441, 4851 |, 14157 |, 2760615 
= 2741-5 - 22 z ef — zi её = E 59 е? – 58 е! – 59 е! 
_ 8690825 ys 112285459 a», TA 
2 2 


证 明 : 椭圆 的 参数 方程 是 X=acos0,y=bsin0,(0<0<27x). 
若 某 条 光滑 曲线 ， 能 用 参数 方程 表示 : x= 六 (?),y=Y(1)，Q <1< p, ZRAKE 


REA: L= VLX OF +[Y O d моти 


L=4 Б Va? sin? 0 +b? cos? 040 


= afè a° (1—cos’ 0) +b’ cos’ 040 


=4|: V1-e cos? 040 
2 p2 
ке, e= | E = 是 离心 率 
a a 
利用 级 数 公 式 : 
1 1 1 
人 (—-n+l) А Н 
Мах=1+у2 2 2 руу SM рех <), 
== n! = (2n)!! 
a 2n-1)!! z 7 ТР 
А f? sin” tdt = Cn DU! 7 ， 借 助 逐 项 积分 法 可 求 出 大 的 无 穷 级 数 表达 式 ， 
0 (m! 2 
L=4a [ 21-е? sin? tdt =27T (a,b) 
(2n- ШЕ е? 
-2лай- 
ла > Onl ro 
lp 141 276061 
s a Naas те-% ее et et e 一 ae 56 
8690825 |, 112285459 , 
22 е 256 е =). 
1. 


C У ъа Š 


引 理 3. \2a?+2b? = 2a[1+ 22 2 сеу 
n=l n. 


l, 24 2. 40 56 w 336 > 528 , 54912 
240—226 6 ве gH 660—506 те аа 6 
91520 , 622336 220 


–..). (16) 


= 232 e 236 


2 
正明 : 2a? +2b? = J2a ІНЕ; =V2a 1-2 
a 


= 
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Е 
= 241+ ——)" 
[ > = ( zA 
1 2 2 40 56 336 528 54912 
05 е =, C =a ЕГІ т: et -ma ere 230 е" 
91520 622336 
22 e = 236 E ==. sy, 
1 1 1 
(二 一 | 一 一 1 十 1 
O ) G ) 


| ce] 
п. 
1,4, 8 ¢ 30, 896 » 10752 ,, 33792 , 7028736 
= 2al -ze Ea E ве" = 58 е!? 58 е!“ zm 59 
23429120 ,。 318636032 ,, 
22 67-0; 6 


ШЕ НЯ: RAE!) =a(1+Vl-e’) 
a 


16 
e 


e). (17) 


P EN Gane) 
l (-еу1 

n: 

4 8 320 896 10752 33792 7028736 

25 e = 28 e= 24 ee 216 т 220 eve 222 = 2% е? 

23429120 , 318636032 5, 

22 G es 256 е-е). 
а зр ТЕР Cael 

引 理 5. 2Vab =2a[l+ 4 4 4 (-e2)"] 
n=l n! 


2-62 14 6 616, 1848 о 23048 ,, 76912 „ _ 16612992 4 


-2а(1- > e — 


= 2а(1 т 2 е 26 25 aie 216 2% 2% 2% 
57222528 801115392 
22 ers 236 te” —..). (18) 
iM: 2Vab =2a, | aie 
a 
1 lÍ 
2 rae jen (nD 
=2all+ >` (—e2)”] 
n=l n. 
1 6 14 616 1848 23048 76912 16612992 
-2а(1- 22 Е = 26 ae 28 en 214 ens 216 Ета 220 a 22 ee 230 е" 
57222528 |, 801115392 „ 
22 Co 256 e =). 
iee n n S, 
a+b mi 2 (n+])! 
1 8 20 896 2688 33792 109824 23429120 
=24(1- 22 E 26 в 28 ao 24 es 216 ar 220 a 222 a 29 e° 
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Е тот os 1 — > ). (19) 
д 
> 2 о 
正明 ， 4ab de l-e 4а! е aa 
a+b 1 b 1+V1-e” e 
a 
i 1-13 ee (2n-1) >] 
m 2 (n+l)! 
1 8 20 896; _ 2688 33792 109824 23429120 
=24(1- 2 се 26 с 28 g= 244 216 2 = 220 ere 222 p= 29 е" 
76959008 |, 1100742656 5, 
= ез- геу: 
222 2% 


3. 主要 结论 与 研究 过 程 
椭圆 周 长 公 式 一 直 是 数学 史上 一 张 缺 了 角 的 纸 , 清 代数 学 家 项 达 明 曾 推 到 十 分 精确 的 项 
达 明 公式 四 ， 


1.3 , LL.32.5 


1 
L=2na(l- k ШОҒЫ Зан > O (20) 
гер 
=== RHA, а. b 为 半 长 轴 和 半 短 轴 . 


BA \ 起 的 一 个 积分 证 明 见 引 理 2 很 显然 ， 这 个 公式 十 分 繁琐 ， 
文 [人 、[8 中 ， 引 入 了 a,2 的 平方 平均 、 算 术 平 均 、 几 何平 均 、 调 和 平均 ， 并 按 离 心率 
的 正 整 次 容 展 开 窜 级 数 ( 见 引 理 3-6) 考虑 到 这 五 个 展开 式 之 中 都 是 形 如 
a(l- xk’ —x,k* – х, ++) (21) 
的 形式 ， 不 妨 将 其 各 项 系数 的 相反 数 归纳 成 下 表 ， 为 了 方便 比较 ， 同 时 附 上 7T(a,D) PAK 
ak” 的 各 项 展开 系数 的 相反 数 〈 见 下 表 ). 


k WJ; T(a,b) S(a,b) A(a,b) G(a,b) H(a,b) 
2 1 1 1 1 
2 2 22 27 22 
4 3 2 4 6 8 
2° 26 05 926 0. 
6 5 2 8 14 20 
pa 85 oe 2% 25 
8 175 40 320 616 896 
pa РЫ РАН са oP 
10 441 56 896 1848 2688 
26 216 216 216 216 
12 4851 336 10752 23048 33792 
220 2% 2% 220 220 
14 14157 528 33792 76912 109824 
222 22 222 222 222 


10 
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16 2760615 54912 7028736 16612992 23429120 
2 — 230 230 230 230 
18 8690825 91520 23429120 57222528 76959008 
2% 292 232 2% 2% 
20 112285459 622336 318636032 801115392 1100742656 
236 236 236 226 236 
TE: 表 中 的 系数 从 大 项 到 了 项 可 在 文 [7 了]、[8] 中 查 到 、K2 项 到 大 2 项 由 maple 软 件 算 


{des 
Co 


АГУ T (a,b) 5 (а,Ь). A(a,b). G(a,b). H(a,b) 的 不 等 关系 .不妨 设 

T(a,b) = mS(a,b)+nA(a,b) + pG(a,b)+qH (a,b). (22) 

考虑 到 精度 要 求 ， 四 元 方程 组 列 出 四 条 方程 ， 则 有 
m+n+pt+q=l, 


2т+4п+6р +87 = 3, 


(23) 
2т+8п+14р + 204 = 5, 
40т + 3200 +616р+ 8964 = 175. 
coon 3 9 5 
Ж т. п. рї, 人 
考 上 处 到 该 方程 组 中 m 、9 项 系数 之 和 等 于 n、p 项 系数 之 和 ， 运 用 maple 软 件 解 得 


3 9 X S R 
m= р" Е L 16 -gg=9.《〈 该 方程 组 有 多 组 解 的 原因 参照 后 文 推广 部 分 ) 


94-0, x 
T(a,b) = a= Sta, b)+? A(a,b)— =. нда, Ь)]. (24) 
借助 (24) 1. are R ~ ee ears 即 
L< 2а 3 568, b)+— = Ala, b)-= Ha b)]. (25) 
— el Xi TM LA 组 强度 十 分 之 高 的 上 下 界 : 

ET 
= (<)L(<) = терца +b)-4Jab -2a (26) 
即 


--1234(а,і) ~5H(a,b)-2S(a,b)|< T(a,b) < 2021464,6) ~2G(a,b)-3H(a,b)]. 


27 
检验 这 组 上 、 下 界 ， 可 以 得 出 他 们 的 系数 ， akraq нанне, ie 
сол, wR НТ 850 
ee ші. тян 
Ta DD) = 2 A(a,b) ~~ H(a,b)-—S(a,b)- = G(a,b). ТР 


Т. 441 1213 1771 
22 26 25 ks 214 k° 216 k” 218 k”? 2/7 ke 


-а(1 
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172843 ,1。 544885 18 1762973 
gë 226 大 一 228 k 230 


开 得 起 有 项 素数 为 , AFT (a,b) 的 展开 项 系数 . 考虑 到 (27) =Q yt k 


A 


k... (29) 


项 系数 头 


， 将 (27) 上 界 与 (28) 二 分 ， 求 出 


7 3 1 
T(a,b “SA ,b) -—H (a,b) -—G(a,b) -—S(a,}). (30) 
(a,b) КаРа НО а СОБ) 2500) 
cen s 
1 3 5 175 441 4851 14157 

(30) = а( 22 k’ 26 3 28 k‘ 24 к 216 k” 220 ke 22 k“ 
=e 6 — рв ae iad ais 

对 比 可 知 , 它 精确 到 “项 , FORE k'° TG L BJ k'° 项 系数 作 商 , 比值 均等 于 0. 99998, 
HEM AY TS Lik’? 项 系数 之 比 ， 亦 高 达 0. 9999. 通过 maple 软 件 进行 验算 ，(30) 式 是 
Т(а,Б) й ЕЯ. 

那么 ， 又 是 否 能 找到 一 个 同样 精确 的 下 界 呢 ? 

假设 推导 得 出 的 下 界 与 7(a,5) 作 差 ， 大 项 系数 小 于 0. Ek D k 项 一 直到 大 无 穷 次 
方 项 系数 均 小 于 0， 那 么 该 下 界 精确 度 不 足 ， 只 需要 另外 寻找 一 K 项 系数 与 7(a,5b) 相同 的 
下 界 即 可 得 到 更 加 精确 的 结果 ;， 若 这 个 差 的 E>” 的 系数 大 于 0， 则 可 能 在 磊 取 某 一 个 值 的 时 
候 差 大 于 0， 这 求 得 的 就 不 是 界 了 . 因此 为 了 提高 下 界 的 精度 ， 需 要 让 下 界 与 7(a,b) 的 差 更 
多 项 4. 0. 

Б tr eR MIE AG Е (230 的 关于 (m,n, p,q) 的 解 集 : 


3 9 5 484 
"160078 00716 00. 4, 将 gq = 0 代入 ， 得 到 k? 项 系数 为 一 


-代入 ， (| kU 项 系数 为 


4852 


assi 
220 , 将 9 = LRA, 得 到 kk 项 


由 此 猜想 : We 的 2 项 的 系数 为 sat 将 整数 i 
ВЕРЕ pss "e Q 时 猜想 也 成 立 . HEB) T (a,b) 展开 


项 当中 k mU жад ， 据 此 有 推论 ， 假 设 实数 组 (m,n, p,q) 使 得 


тб(а,Б)-пА(а,5)-рСе(а,Б)-аН(а,Б)> ал Т(а,5), ЯЯАДЫЕНҒ4 = БІ . 


32 
、 ik 4 14 14168 16 ур. ` 
进一步 考察 (29), (30), (29) А 项 的 系数 为 一 28? k ° 项 的 系数 为 
— 2709488 ，(30) 的 15 项 的 系数 是 22760560 
2 


2765488 26065. 29090. UN (29), (зо) 式 连 续 多 次 进行 二 分 


2% 230 230 
2760615 1 
法 ， 有 其 中 菜 次 结果 的 K 项 的 系数 小 于 -一 50- ， 而 下 一 次 二 分 的 结果 大 于 


2703488 
т 259 Š 
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经 计算 ， 利 用 (29)、(30) SRE 6 次 二 分 ， 即 g = 一 时 ， 有 一 个 十 分 精确 的 
Pie AGE Ge eC A (32) 
2048 2048 2048 2048 
泰勒 展开 ， 有 
2753 449 65 191 
—— А(а, Б) -—— H(a,b) -—— G(a,b) -—— S(a,b 
2048 (4,0) 2048 (0) 7048 0 ) 2048 GP) 
1 3 5 175 441 310465 906059 2760637 
=] pe pF же rod r > 12 _ 5 14 _ 2% kê... 
же Ey ONP м^ 897 S| 48 = Wee 、 
同 理 第 七 次 二 分 ， de 可 以 得 到 男 一 个 同样 精确 的 近似 值 : 
920952 р P ады CG De Sa (33) 
4096 4096 4096 4096 
泰勒 展开 ， 有 
5505 897 383 129 
— A(a,b) -—— H(a,b) -—— G(a,b) -—— S(a,b 
4096 (2,0) 4096 ap) 4096 ) 4096 (59 
1 3 5 175 44] 620929 1812107 
=1 2 k? 26 k 25 k° 24 к 216 к" 27 ke 25 ke 
5521197,1; 2172555 ,1。 898081479 | 5, 
站 k sm k т K... 
2 2 2 


至 


接 下 来 ， 考 察 (28), (32), (33). 
H (28) 及 引 理 1、2, 算出 两 个 椭圆 C :a=1,b=0.05, С,:а=1,6= 0.5, WAK, 


并 与 文 [10] 当 中 椭圆 周 长 的 真 值 作 比 较 . 可 以 发 现 ， 由 (28) 算出 C, 的 周 长 大 于 文 [11] 中 提 
I 的 真 值 , 但 算出 CC; 的 周 长 却 小 于 文 [11] 中 提 到 的 真 值 . 这 也 说 明了 (28) 在 某 些 时 候 大 于 
7T(a,b)， 某 些 时 候 小 于 7(a,5b)， 故 (28) 并 不 是 7(a,b) 的 上 界 或 下 界 . 同 理 考察 (33) 可 


以 得 到 一 样 的 结果 ， 故 (33) BRET (a,b) 的 上 界 或 下 界 . 


再 考察 (32), (32) 泰勒 展开 式 中 ，K”、k 项 的 系数 均 小 于 7T(a,b)， 但 是 (32) = 


的 "项 的 系数 却 大 于 7(a,5)， 这 说 明 (32) 在 某 些 情况 下 大 于 7(a,5) ， 某 些 情况 下 小 于 


H (28) 以 及 (32) 知 ， 仅 仅 由 泰勒 展开 式 的 前 面 几 项 的 系数 不 能 求 出 下 界 . 为 了 求 出 下 界 ， 


不 妨 设 a =1,b=XxX， 则 将 近似 公式 化 为 


3 l+x 9 l+x 5 2x 
----4)х + x + x Vx+ x А 
Се 4) б G q) 5 (тє 4)хУх+9 TF 
3 l+x 9 l+x 5 2x 
令 = (一 十 十 十 十 
Ә-С Әм + G- Ot q)x Vx qx 
ea) x +1 5 
一 == 
' 16 2 9 16 2q 
则 有 f (x)= + + + . 
f) yl aa Guy 


& f (x) =0, TAM maple 解 得 其 中 的 实数 根 x 为 和 gq 有 关 的 非 负 实 根 ，9 取 不 同 值 


的 时 候 x 也 有 不 同 的 值 . ZE q = О 的 时 候 有 x ғ 0.01920312361. 


由 此 可 以 说 明 , 当 0<xz<x 时 , MB f(x) <0, W f(0) > f(x). FBT, x) fE 


x €[0,1] 上 显然 是 单调 递增 的 ， 有 7T(,0) < T(x) .又 因为 在 0<x<1 时 f(x) <T(,x)， 


所 以 0<xz<w ЕНЕЯ /(х)<Т(,х). ik f(0) < T(1,0) .在 此 的 基础 上 ， 使 用 maple 软件 进 
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行 更 加 精细 的 计算 ， 可 以 得 出 7(c,DD) h 
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的 下 界 : 


BOG ED SG (34) 

2048 2048 2048 2048 

考虑 到 运用 二 分 法 可 能 会 求 出 满足 f(0) = T (1,0) 的 近似 公式 ， 本 文 不 在 (34) Я (30) 

的 基础 上 继续 采用 二 分 法 . 

fe PR, BAI GO AM (34) 式 进 行 误差 分 析 : 

а |b 下 界 下 界 / 真 值 “ 真 值 ” 上 界 / 真 值 LH 
(34) (30) 

1 | 0.00 | 3. 999438207 | 0. 9998595518 | 4 1. 020668759 | 4. 082675036 
1 | 0.05 | 4. 009009874 | 0. 9974086482 | 4. 019425620 | 1. 002876644 | 4. 030962381 
1 | 0.10 | 4.057412652 | 0. 9983854405 | 4. 063974180 | 1. 000999318 | 4. 067961130 
1 | 0.15 | 4. 122198813 | 0. 9990544382 | 4. 126100296 | 1. 000378119 | 4. 127660157 
1 | 0.20 | 4. 199743250 | 0. 9994608155 | 4. 202008908 | 1. 000152428 | 4. 202649005 
1 | 0.25 | 4. 287922836 | 0. 9996996995 | 4. 289210888 | 1. 000062318 | 4. 289477794 
1 | 0.30 | 4. 385194345 | 0. 9998368131 | 4. 385910068 | 1. 000025327 | 4. 386021082 
1 | 0.35 | 4. 490352261 | 0. 9999137534 | 4. 490739572 | 1. 000010111 | 4. 490785007 
1 | 0.40 | 4. 602419414 | 0. 9999558717 | 4. 602622520 | 1. 000003928 | 4. 602640601 
1 | 0.45 | 4. 720586687 | 0. 9999782750 | 4. 720689244 | 1. 000001465 | 4. 720696159 
1 | 0.50 | 4. 844174652 | 0. 9999897903 | 4. 8442241102 | 1. 000000518 | 4. 844226618 
1 | 0.55 | 4. 972606862 | 0. 9999954675 | 4. 9726294004 | 1. 000000171 | 4. 972630249 
1 | 0.60 | 5. 105390207 | 0. 9999983222 | 5. 1053997727 | 1. 000000051 | 5. 105400035 
1 | 0.65 | 5. 242099893 | 0. 9999992929 | 5. 2421035996 | 1. 000000014 | 5. 242103672 
1 | 0.70 | 5. 382367708 | 0. 9999997634 | 5. 3823689815 | 1. 000000003 | 5. 382368998 
1 | 0.75 | 5. 525872670 | 0. 9999999933 | 5. 5258730402 | 1. 000000001 | 5. 525873043 
1 | 0.80 | 5. 672333493 | 0. 9999999851 | 5. 6723335778 | 1. 000000001 | 5. 672333578 
1 | 0.85 | 5. 821502467 | 0. 9999999977 | 5. 8215024803 | 1. 000000001 | 5. 821502480 
1 | 0.90 | 5. 973160432 | 0. 9999999999 | 5. 9731604325 | 1. 000000001 | 5. 973160433 
1 | 0.95 | 6. 127112637 | 0. 9999999999 | 6. 1271126366 | 1. 000000001 | 6. 127112637 
1 | 1.00 | 6. 283185307 | 0. 9999999999 | 6. 2831853072 | 1. 000000001 | 6. 283185307 

YE 1: 这 个 表 的 上 下 界 由 于 是 学 生计 算 器 计算 而 得 ， 只 能 取 10 位 数字 ， 其 中 最 后 一 位 


是 四 舍 五 入 求 得 故 上 界 在 5 接近 a 时 可 能 由 于 四 含 五 入 而 导致 显示 的 上 界 的 数值 比 “ 真 值 
略 小 ， 而 下 界 的 数字 比 “ 真 值 ” 略 大 ， 这 是 属于 学 生计 算 器 不 足 而 导致 的 ， 并 非 上 界 与 下 界 
不 精确 而 导致 的 . 
注 2: 表 中 的 “ 真 值 ” 见 文 [10]， 文 [10] 中 不 精确 部 分 由 maple 软件 算出 . 
考察 上 表 可 以 得 出 : 当天 越 小 时 ， 上 下 界 的 精确 度 越 高 . 从 几何 角度 来 说 ， 离 心率 越 小 ， 
椭圆 越 接近 圆 ， 而 本 文 的 上 下 界 是 将 椭圆 类 比 成 圆 推导 出 来 的 ， 因 此 椭圆 越 接近 圆 ， 上 下 界 
WA BOR UL ASAE. 
wa=1, b2J[0.1]! 
到 以 下 几 个 结论 ; 
I 、 对 于 所 有 椭圆 ， 用 2x x (30〉 作 为 近似 公式 计算 
Ci) 90% 以 上 的 概率 使 得 其 误差 小 于 干 分 之 一 . 
Cii) 75% 以 上 的 概率 ， 使 得 其 误差 小 于 万 分 之 一 . 
(111) 略 低 于 65% 的 概率 ， 使 得 其 误差 小 于 十 万 分 之 一 . 
Civ) 50% 以 上 的 概率 ， 使 得 其 误差 小 于 百 万 分 之 一 . 
(у) 40% 以 上 的 概率 ， 使 得 其 误差 小 于 于 万 分 之 一 . 
H. Ot PRA, НІ2лх (34) 作为 近似 公式 计算 
Ci) 85% 以 上 的 概率 使 得 其 误差 小 于 干 分 之 一 . 
Cii) 超过 65% 的 概率 ， 使 得 其 误差 小 于 万 分 之 一 . 


的 某 一 随机 值 ， 对 上 表 分 析 (30) 式 和 (34) 式 的 误差 ， 可 以 得 


HAK, A 


其 周 长 ， 有 
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Gii) 超过 45% 的 概率 ， 使 得 其 误差 小 于 十 万 分 之 一 

Civ) 35% 以 上 的 概率 ， 使 得 其 误差 小 于 百 万 分 之 一 . 

Ш. (30) RA GO RH T(a,b) 的 误差 概率 同 工 、 工 ， 

考虑 到 (30) 式 的 泰勒 展开 式 可 以 精确 到 大 的 14 Tae, Th (34) 式 的 泰勒 展开 式 仅 
逼近 但 不 精确 到 左 的 12 WIT, (30) 会 比 (34) 略微 精确 . 

回 到 Toader 平均 ， 则 有 E 


2883 579 195 
2048 A(a, aie (a, ШЕСІ (а UET b) <Т(а,5) 
< Aab) -Z На, b)-= Gla, b)- 二 St 5. (35) 


查阅 参考 文献 ， кан EE) 的 趟 等 : 式 有 许多 . 
接 下 来 ， 本 文 对 比 上 文 发 现 的 上 、 下 界 与 参考 文献 中 提 及 的 上 、 下 界 . 
ГАТ ЖТТ (а, Б) 的 不 等 式 


= 2 = 2 
ДЕТЕ ал кі Е aie eta 
2 I+ 1-2 8 
2 2 
人 (36) 
л 1+1? л 
ма? -b 


中 
a 
BSCE, (36) 式 与 (5) 式 等 价 . 
用 定理 1 中 得 出 的 上 界 与 5) 式 中 上 界 作 差 ， 利 用 maple 软件 计算 得 出 的 结果 如 下 : 


4337-32 0% ke 2-2. 89 2 
8 1+41-k? +1-k° 6 
-da[( 一 - +(2-—)-(l¢+V1-k? 
a 2): Т БЕГІН аа +( a (+ )] 
2 = by 


- а, b)- Hla, b)-= Gta, b)- Sa »b)— ib 
a 


- 4, b)- Hla, b)-= Gta, b)- S(a, 5)- [= A(a, b)+0-Ż)H(a, b)] 


4 4 
P -Sra b)-= Gla, b)- S(a, b) 
з ` a D. 1. zi Еее 

4 X 4 A 4r 64 256 8 64 л 

465 _ 5 1 At 7 tz 7 ,1351 zn 

16384 647 64 л 512 л 1287 65536 
140.16653 ，3 Ate 21 a (53427 165 Atr 33 yu 
1048576 256 л 512л 4194304 16384 л 10247 


(706057 2 Arr 49 ae 
67108864 16384 л 16384r 268435456 131072 л 327687 
32813523 221 -4+л 2431 py 

4294967296 32768 л 1310727 


2385515 1001 -4-л 715 ,ig 
+ k 
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114464493 u 12597 —4+z 4199 ТЕ 
17179869184 2097152 ~x 2621447 
ц 808069563 u 11305 -4ғл 29393 Jk” 
137438953472 2097152 ~x 20971527 


2880773483 P 81719 -4+л 52003 


+( у 
549755813888 16777216 л 4194047 
41429163149 37145 —4+m 185725 ТЕ 
8796093022208 8388608 л 167772167 

” 150042970683 4345965 —4+m _ 334305 ТЕ 


35184372088832 1073741824 л 335544327 

+ 4375026633141 N 3991995 -4-л 9694845 ТЫ 
1125899906842624 1073741824 л 10737418247 
16035710388339 29464725 -4-л 17678835 p“ 


二 一 一 一 一 一 一 一 十 一 一 一 一 一 一 一 一 一 一 一 
4503599627370496 8589934592 л 21474836487 
қ 236271807547479 А 3411705 -4-л 64822395 ТЫ 
72057594037927936 1073741824 л 85899345927 


405992895 -4-л 119409675 i 874178618994681 


137438953472 ~x 171798691847 288230376151711744 
6494516919141405 y 378699255 -4-л 883631595 


2305843009213693952 137438953472 xm 1374389534727 


+( Jk” 


+( Jk +...] 


(37) 
J maple 软件 将 (37) RAF ЕАС Pra BSR, RM (37) 式 项 直 
2. , Mk 项 开始 系数 为 正 值 . 又 因为 玉兰 0， 所 以 定理 1 中 的 的 上 界 与 
(5) 式 不 分 强 弱 . ІН maplel6 软件 更 加 精细 地 计算 得 : 


当 0 << 0.0399212204… 时 ， 


et 


Ala 1b) Ha, b)- Gla, b) 5 Sta be- 2 aab+(a-by ， 
a+b 
当 » >0.0399212204---14, 
а 
Z ада, b)- Hla, b)-= Gla, 已 -本 5 S(a pes 2 D. 
a 


这 说 明 ， "4 > 0.0399212204.… 时 ， (5) AE GD 式 的 上 界 弱 . 假设 a =1, b » [0,1] 
a 


中 的 某 一 随机 值 ， 可 以 得 出 结论 : (35) 式 的 上 界 比 (5) 式 的 上 界 更 强 的 概率 约 为 96%. 也 
就 是 说 ， 在 绝 大 部 分 情形 下 ，(35) 式 上 界 比 (5) 式 上 界 更 强 . 
同时 ， 用 定理 1 中 的 下 界 与 (5) OPA RAE, ЖІН maple 软件 可 以 得 


412583. l+Vl-k? 289 2(41- k’ -1+k°) 6l AE 195 ІШ a 


2048 2 2048 k? oer 2048. 
_ 2 < 2 
oo К +1-k а A-k] 
2 1-41-42 8 
2883 579 61 195 5а? + 6ab + 5b’ 
=——— A(a,b) -—— H (a,b) -—~ G(a, b) - —_ S(a,b) -—___— 


2048 2048 2048 2048 8(a+b) 
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2883 579 61 195 5 1 
Se 

GON аттты Бов ol K s Р 
= ES ЖЕЛЕ 71) БЕЛЕ) ШЕ а 

2048 2048 2048 2048 
СЫ 7 ооф 8893 іш, 3867 и, 157993 |, 


16384 65536 67108864 268435456 1073741824 
78001 is, 38710377 | 148541601,» 
536870912 274877906944 1099511627776 
566769259 ы, 539077519 1,20) e5 
4398046511104 4398046511104 
考察 每 一 项 的 系数 , WAT 0, ХК 20, 所 以 定理 1 中 的 下 界 亦 更 加 逼近 7(a,D) . 
因此 ， 我 们 可 以 得 到 一 串 不 等 式 甸 9: 


КЕ 2 _ 2 
ек о pra 
2 l+V1-k? 8 
Шер A он) 42664 Sun 
2048 2048 2048 2048 
<T(a,b) 
43 7 3 l 
22 A(a,b)-— H(a,b) —— G(a,b) —— S(a,b 
ЖЕ (4,0) — 35 GB) 3 9) 55 (0D) a 
< 39 
2 2 
E a a +1-k +2 DP ha) 
л 1+V1 一 大 л 
212 
раа р 
а 
与 此 同时 ， 文 [9] 中 也 提 到 
“sr Ел), (40) 


A 
同样 的 ， 与 (5) 式 作 比 较 ， 可 以 得 至 
a+b К 5а? + 6ab +5b° 


| 


< <Т(а,5). (41) 
2 8(a +b) 
НУ. 5a’ +6ab+5b° a+b = 5a’ +6ab+5b° —4(a+b)’ = (a—by a 
8(a +b) 2 8(a+b) 8(a +b) 
ia mab +(a—b) 2 4а-ль 不 分 强 弱 ， 
A a+b 2л 
жасыр 4азлЬь 2 mab+(a-by _(4a+zb)Xa+b)-4[zab+(a-b)] 
2277 OR т. a+b g 2z(a+b) 
_ (12—3x)ab+(x —4)b? _ b(á —z)(3a —b) 
g 2л(а+Ь)  2m(a+b ` 
“43a > b W, 4а-ль 、2 zab+(a-by) | 


л л a+b 
4a+zb 2 _mab+(a-by l 
2л A a+b | 


当 3a = Ьі, 
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ДЕ 2 
当 34 < b itt, 4a+zb 2 zab+(a b) l 
2л A a+b 
H L =2zT(a,b) Ka > b 可 得 到 
2 2 ру 
sara e +6ab+5b 12у < Azab + (a b) atar 
4(a+b) a+b 
不 难看 出 Ш» 定理 1 与 定理 2， 确实 很 强 . 
综 上 所 述 ， 我 们 有 一 串 不 等 式 链 : 
2 2 
a+b 1 l+vȚvl-k +1l-k — 
2 2 I+ 1-2 
2883, 579 195 
< --- 一 一 一 Б) -—— S(a,b 
204840227 2048702227 5048 022) 2048500) 
<T(a,b) 
43 


ag b)- SH, b)- = Gla, b)- 5504, Б) 


2, 2 
8 1+1=4 +1-k арусе 
Ж т 


I (= 2)x +(2 


z 1+V1 一 大 
相应 的 椭圆 周 长 公 式 有 
1+y1-k° +1—- 2 
m(a+b) < zal + + k°)] 
l+V1-k? 
<j 28834(a, Б) -579H(a,b)—61G(a,b)—195S(a,b)] < L 
T [434(a,b) -7H (a,b) =3G(a, b) = S(a,b)] 
ЕТА sS k tlk’ +(47 -6)1+ 31-2) <4a+ zb 
1+ 1-2 


(42) 


(43) 


(44) 


推广 : 在 研究 T(a,5) 近似 值 的 过 程 当中 ， 发 现 了 很 多 组 近似 公式 的 泰勒 展开 项 系数 一 


直到 左 的 10 W ABE T (a,b) 的 泰勒 展开 项 系数 相同 . 但 是 用 待定 系数 法 列 出 四 元 一 次 方 
程 组 并 解 方程 的 角度 来 说 ， 这 个 四 元 一 次 方程 组 应 该 具有 一 组 解 . 这 个 结果 明显 与 研究 过 程 


的 情形 不 同 ， 这 是 为 什么 呢 ? 在 整个 研究 过 程 当中 ， 发 现 了 很 多 组 (mm. n. p. q) 使 得 


mS(a,b)+nA(a,b)+ pG(a,b)+ gH (a,b) ВЕ ЖИІ AY 10 KE, (ASS AIL RA 
7 (т. m p. а) ИР 14 mS(a,b)+nA(a,b)+ pG(a,b) + gH (a,b) 精确 到 大 的 14 
Хә» BEE ESC ACE A SE ee BJ PRE Ae IPL, 大 的 12 RERA T (a,b) Pk 


的 12 W: At AN AP. 
m+n+ p+q=1 
2т+4п+6р +87 = 3 


2т+8п+14р + 204 = 5 
40m + 3200 + 616р + 8964 = 175 


这 个 方程 显然 有 解 . 同时 有 : 在 同一 条 式 子 中 ，m 的 系数 与 4 的 系数 之 和 等 于 n 的 系数 


与 p 的 系数 之 和 .不妨 设 (m. п. p. q) 为 该 方程 组 的 世 


其 I 


一 组 解 . < 


m, =m, + Ax,n, =n, — Ax, p, = p, —Ax,q, = 4 + Ax(Ax + 0) 代入 原 方程 组 ， 显 然 成 立 ， 
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ы 


МЕ (т, п, Py» 4) 也 为 该 方程 组 的 解 . 由 此 可 以 证 明 , 这 个 方程 组 有 无 数组 实数 解 . 
am+bn+cp+dq=e, 
a,m+b,n+c,p +d,q =e, 


此 可 以 大 胆 猜 测 : 已 知 方程 组 有 实数 解 . 


a,m+b,n+c,pt+d,q = e, 
a,m+b,n+c,p +d,q =e, 
HiX i =1,2,3,4 1а, +b, =c, +d, WZI FERA ABA ЭСТЕ. 
更 一 般 的 ， 有 以 下 结论 : 由 实 条 两 两 不 同 的 方程 组 成 的 2N 元 一 次 方程 组 中 ， 若 其 中 
个 未 知 数 在 任 一 条 方程 中 的 系数 和 均等 于 另外 和 个 未 知 数 在 这 条 方程 中 的 系数 和 , 那么 该 方 
程 组 有 无 数组 实数 解 或 者 无 实数 解 . 


4. 椭圆 周 长 近 似 公式 的 应 用 

椭圆 积分 、Toader 平均 ， 以 及 椭圆 周 长 近 似 公 式 在 数学 、 物 理 、 天 文学 以 及 工程 中 有 
着 广泛 应 用 .下面 从 文 [11] 中 的 几 道 习题 略 谈 定 理 2 的 应 用 . 
例 1.( 选 修 2-1P49 页 习题 2-2A 组 第 9 题 ) 彗星 “紫金 山 一 号 ?是 南京 天 文 台 发 现 的 . Е 
的 运行 轨道 是 以 太阳 中 心 为 一 个 焦点 的 椭圆 . 测 得 轨道 的 近日 点 (距离 太阳 最 近 的 点 ) 距 太阳 
中 心 1. 486 天 文 单位 ， 远 日 点 (距离 太阳 最 远 的 点 ) 距 太阳 中 心 5. 563 天 文 单位 (1 天 文 
单位 是 太阳 到 地 球 的 平均 距离 , 11. 5x10km)， 且 近日 点 、 远 日 点 及 太阳 中 心 在 同一 条 直 
线 上 ， 求 轨道 的 方程 . 

解 : 由 天 文学 知识 知 , 太阳 直径 1, 392, 000 公 里 (地 球 直径 的 109 倍 ) 可 以 换算 成 0. 0928 
天 文 单位 . 据 此 可 以 求 出 该 彗星 运动 轨迹 的 长 轴 为 3. 571 天 文 单位 , 短 轴 为 2. 899 天 文 单位 . 

2 2 
不 难 求 的 轨道 的 方程 是 一 二 + 一 =1，( a,b 的 单位 均 为 1 天 文 单位 ) 
3.5717 2.899 

由 此 ,我 们 借助 定理 2, 得 到 该 人 造 地 球 卫 星 的 轨道 周 长 为 20. 38095945 至 20. 38095966 
天 文 单位 , 可 见 定理 2 的 精确 度 很 高 . 

B] 2. (选修 2-1P49 页 习题 2-2A 组 第 10 Ж) 已 知 地 球 运 行 轨道 是 长 半 轴 长 
а=1.50х10 km， 离 心率 e = 0.0192 的 椭圆 ， 且 太阳 在 这 个 椭圆 的 一 个 焦点 上 ， 求 地 球 到 太 
阳 的 最 大 和 最 小 距离 . 

Ж. Ша-1.50х10%ш, 得 b=Va? -c =avl—e? =1. 4997х 10%, 

由 此 ， 我 们 借助 定理 2， 得 到 该 地 球 轨道 的 周 长 为 

9. 423835509 х 10°km< (==) L< (=) 9. 423835509km. 可 见 定理 2 精确 度 很 高 . 

例 3. (选修 2-1 复习 参考 题 P80A 组 第 1 题 ) 如 图 , 我国 
发 射 的 第 一 颗 人 造 地球 卫 星 的 运行 轨道 , 是 以 地 心 (地 球 的 中 
心 )F2 为 一 个 焦点 的 椭圆 , 已 知 它 的 近地点 A( 离 地 面 最 近 的 点 ) 
距 地 面 439km, 远地点 B 距 地 面 2384km. ЗЕН. F2, А. В 在 同一 
HRE, 地球 半径 约 为 6371km, 求 卫星 运行 的 轨道 方程 (精确 
到 1km). 

Ж. 以 直线 AB 为 x 轴 , 线段 AB 的 中 垂 线 为 了 轴 建 立 

如 图 所 示 的 直角 坐标 系 ，AB 与 地 球 交 与 C. D 两 点 . 


ARMIA 207—1 (a>b>0), 

由 题 意 知 : |AC|=439, | BD |=2384, 

ЕСЫ E,D 6371. 
WU): a-c = |04|-|ОЕ,|=|Б,А| = 6371+ 439 = 6810, 
a+c=|OB|+|OF,|=|F,B| = 6371+ 2384 = 8755. 
解 得 a = 7782.5,с = 972.5, 2b 7722. 
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ba 


х? 


2, 
故 卫 星 的 轨道 方程 为 约 一 十 一 二 一 1. 
77835 7722 


由 此 ， 借 助 定理 2 计算 得 到 该 人 造 地 球 卫 星 的 轨道 周 长 世 满足 
48710.58258km = (<)L(<) = 48710.58258km , НІН L = 48710.58258km. пу Mae HE 2 
的 精确 度 很 高 . 

有 兴趣 的 读者 ,还 可 以 借助 定理 2 研究 人 类 已 发 现 的 各 种 下 星 ( 见 文 [12]) 的 轨道 周 长 . 


5. 课题 研究 展望 
本 文 提供 了 一 种 探索 Toader 平均 的 方法 ， 也 给 出 了 一 对 十 分 简洁 而 又 较为 精确 的 上 下 
Ӯ, 也 就 是 


2883 579 195 

~ АРЫ т. ea 2 
2048 422) 30487020) On K OA 
<T(a,b) 


<5 да у На, b)-= Gta, b)- Sta, b) 


айша, E aN on 调和 平均 来 逼近 Toader 平均 的 方式 
或 许 能 够 应 用 在 求 其 他 不 能 用 初等 函数 表示 的 平均 之 中 ， 甚 至 也 可 以 扩展 至 与 
а,а,,а,, 4, 有关 的 更 高 阶 的 Toader FI 7 (а,,а,,а,,...,а,). 也 可 以 对 a、b 进行 加 权 ， 
求 得 其 四 个 加 权 平 均 数 ， 再 用 同样 的 方法 来 通 求 7(a,D) . 

本 文采 用 的 四 种 平均 值 逼 求 的 7(a,D) 的 方法 仍 存在 一 些 不 足 ， 如 上 文 所 提 到 的 ， 某 些 
公式 的 泰勒 展开 项 系数 比 7(a,p) 小 ， 但 是 在 其 中 一 个 变量 取 到 较 小 值 的 时 候 却 大 于 
7T(a,b)， 由 此 导致 了 下 界 的 强度 不 如 上 界 , 这 个 问题 是 这 种 方法 缺陷 的 一 部 分 . 但 是 瑕 不 掩 
瑜 ， 本 文 得 到 的 简洁 的 上 下 界 依旧 有 极 高 的 价值 ， 由 此 可 以 衍生 出 一 系列 的 不 等 式 , 可 以 与 
伽 马 函数 等 高 等 函数 相 结 合 ， 可 以 求 出 椭圆 的 周 长 并 结合 开 普 勒 天 文 定律 求 天 体 轨 迹 周 
长 …… 本 文 的 成 果 还 有 诸多 应 用 ， 限 于 篇 幅 ， 不 再 袭 述 . 
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